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Abstract 

For a two-dimensional lattice A with n vertices, the number of spanning trees Nst{^) grows 
asymptotically as exp(nzA) in the thermodynamic limit. We present exact integral expression 
and numerical value for the asymptotic growth constant for spanning trees on various two- 
dimensional lattices with more than one type of vertex given in [l]. An exact closed- form expression 
for the asymptotic growth constant is derived for net 14, and the asymptotic growth constants of 
net 27 and the triangle lattice have the simple relation 227 = (2^iri + lii4)/4. Some integral identities 
are also obtained. 
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I. INTRODUCTION 



The enumeration of the number of spanning trees Nst{G) on the graph G was first 



considered by Kirchhoff in the analysis of electric circuits [2|] . It is a problem of fundamental 



interest in mathematics 



n 



n 

k 



and physics [7|, ISj. The number of spanning trees is clo sely 
related to the partition function of the g-state Potts model in statistical mechanics [q, [iO |. 
There are several ways to calculate Nst{G), including as a determinant of the Laplacian 
matrix of G and as a special case of the Tutte polynomial of G jsl- Some recent studies on 



the enumeration of spanning trees and the ca^ 



culation of their asymptotic growth constants 



In this paper we shall present exact 



on regular lattices were carried out in [11 
integrals for the asymptotic growth constant for spanning trees on several two-dimensional 
lattices given in jj]. It was shown in IJ] that integral identities can be obtained with 
different choice of unit cells, where, in most cases, the sizes of the unit cells are different. 
Here we shall show that an integral identity can be obtained with the same unit cell but 
different choice of directions in the calculation. 



II. BACKGROUND AND METHOD 

We briefly recall some definitions and background on spanning trees and the calculation 

nn 

method that we use [3|, ll5[. Let G = (V, E) denote a connected graph (without loops) with 
vertex (site) and edge (bond) sets V and E. We will only consider simple graphs without 
multiple edges. Let n = v{G) = \V\ be the number of vertices and e{G) = \E\ the number 
of edges in G. A spanning subgraph G' is a subgraph of G with v{G') = |V^|, and a tree is 
a connected graph with no circuits. It follows that a spanning tree is a spanning subgraph 
of G that is a tree and hence e{G') = n — 1. The degree or coordination number ki of a 
vertex Vi ^ V is the number of edges attached to it. A fc-regular graph is a graph with 
the property that each of its vertices has the same degree k. Two vertices are adjacent if 
they are connected by an edge in E. In general, one can associate a.n edge weight Xij to 



each edge connecting adjacent vertices Vi and vj (see, for example ll|]). For simplicity, all 
edge weights are set to one throughout this paper. The adjacency matrix A{G) of G is the 
n X n matrix with elements A{G)ij = 1 if f j and vj are adjacent and zero otherwise. The 
Laplacian matrix Q{G) is the nxn matrix with element Q{G)ij = ki6ij — A{G)ij. One of the 



2 



eigenvalues of Q{G) is always zero; let us denote the rest as A(G)i, 1 < i < n — 1. A basic 
theorem is that NsriG) = (l/n) YYi=i -^(G*)! [3]. For a ci-dimensional lattice A with d > 2 
in the thermodynamic limit, A^5'j'(A) grows exponentially with n as n — » oo; that is, there 
exists a constant z\ such that Nst{-^)^ exp{nzA) as n ^ oo. The constant describing this 
exponential growth is thus given by 5|, |6| 



z/^ = lim — In 

n—*oo n 



NsAA) , (2.1) 



where A, when used as a subscript in this manner, implicitly refers to the thermodynamic 
limit of the lattice A. 

A regular (i-dimensional lattice is comprised of repeated unit cells, each containing u 
vertices. Define a{n, h') as the u x u matrix describing the adjacency of the vertices of the 
unit cells n and n', the elements of which are given by a{n, h')ij = 1 if f j G n is adjacent to 
Vj G n' and zero otherwise. Although the number of spanning trees A^st(A) depends on the 



boundary conditions imposed as shown in 



the asymptotic growth constant z/^ is not 



sensitive to them. For simplicity, let us consider a given lattice having periodic boundary 
conditions. Using the resultant translational symmetry for the spanning trees, we have 
a{n,n') = a{n — n'), and we can therefore write a{n) = a{ni, ■ ■ ■ ,nrf) for a (i-dimensional 
lattice. Generalizing the method derived in for lattices which are not /c-regular, Nst{-^) 
and z\ can be calculated in terms of a matrix Ma, which is determined by these a(n), defined 
as 

Ma(^i, ...,9,) = M'j,-J2 a(^)e^"® , (2.2) 

h 

where Mi is the diagonal matrix whose diagonal elements are the degrees h of the vertices 
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in the unit cell and 9 stands for the d-dimensional vector {9i, ■ ■ ■ ,9d)- Then [5|, Il2l | 

^A = -/ M^a(^i,---,^<^)] , (2.3) 

^ J-^ i=l ^ 

where Dj>^{6i, ■ ■ ■ ,6d) = det(MA(6'i, ■ ■ ■ , 6*^)) is the determinant of the matrix Ma- Notice 
that the calculation is not sensitive to the order of the vertex labeling and the choice of the 
directions 6j, I < j < d. 

It is well known that there are only three uniform tilings of the plane by using one type 
of regular polygon in which all vertices are equivalent, or three regular tessellations, namely, 
the square, triangular and honeycomb lattices. If one allows more than one kind of regular 



polygons and still requests that all vertices are equivalent, there are eight more lattices, or 



semi-regu 
fc-regular 



ar tessellations. These are altogether eleven Archimedean lattices which are all 



Iq . If the restriction that all vertices are equivalent is released, an infinite number 



of tessellations is possible, even with just two types of regular polygons, not to mention if 
non-regular polygons are allowed. Even though mathematically it is not possible to cover 
the plane if regular pentagons or heptagons should be presenst, certain arrangements of 
atoms involving irregular polygons, including pentagons or heptagons, do occur in real- 
world alloys and inorganic crystals. In various common tessellations (including the 
eleven Archimedean lattices), denoted as nets, and their occurrences were given. 

For a lattice A which is not /c-regular, it is convenient to introduce an effective coordina- 
tion number ka, defined as the average number of edges per vertex, 

2e(A) 



lim 



(2.4) 



n(A)^oo ?7,(A) 

For a fc-regular lattice, k = k. Furthermore, we know that the number of spanning trees is 
the same for a planar graph G and its dual G*, and the number of the vertices of G* is given 
by the Euler relation v{G*) = e(G) —n + 1. It follows that the asymptotic growth constants 
of G and G* satisfy the relation Q, Q 



Zg* 



(2.5) 



k/2 - 1 

For a fc-regular graph Gk, a general upper bound for the asymptotic growth constant is 
ZGk < In/c 17 1 . A stronger upper bound for Gk with k >3 was derived in |l8, 1^ that 
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NsAGk) < 



nk In k 



where 



(k-l) 



k-i 



By Eq. ([211]), this then yields [1 



[k{k-2)Y. 



zgu < Hh 



(2.6) 



(2.7) 



(2.8) 



III. ASYMPTOTIC GROWTH CONSTANTS 



The asymptotic growth constants 



sidered by several authors 



QBE, 



or the eleven Archimedean lattices have been con- 
While the relation Zhc = Ztri / 2 for the honeycomb 



and triangular lattices is easy to understand due to the duality [12] (cf. Eq. fl2.5p ). it is 
non-trivial to have the relations Zkag = (^tri + ln6)/3 for the Kagome (equivalently (3.6.3.6)) 



lattice and -2(3.12.12) = {ztri + ln(15))/6 for the (3.12.12) lattice given in Our main 

purpose is to calculate the asymptotic growth constants for other common two-dimensional 
lattices where more than one type of vertex occur. Following the denotation given in [l|, 
we shall quote them as net 12 to net 27, where net 12 to net 17 are tessellations with two 
or three regular polygons (including triangle, square or hexagon), while net 18 to net 27 are 
tessellations containing pentagons, heptagons or enneagons. In addition, [l| mentions the 
B net of YCrB4, the B net of Y2LnB6 and a net with 5^.8 and 5.8^ vertices. Let us denote 
;hem as nets 28, 29 and 30, respectively. The figures of these nets are referred to those in 
l|, and the unit cells chosen for the calculation are shown in Figs. [T]131 Notice that the 
polygons are regular or not does not affect the number of spanning trees, so the unit cells 
can be deformed from those in pL|]. In addition to the numerical values of the asymptotic 
growth constants, we shall derive an exact closed-form expression for Z14, and the relation 
^27 = [ztri -|- ln4)/4. A few integral identities will be given by choosing different unit cells 
or directions. We will use the shorthand notations a = e*^^ and f3 = e*^^ for the elements of 
the matrix Ma when the matrix is large. 

A. Nets with regular polygons 

In this subsection ,we consider important nets in which only regular polygons (including 
triangle, square or hexagon) occur with more than one type of vertex. 

1. Net 12 

Net 12 is the combinations of 3^.4.3.4 and 3^.4^ vertices. A primitive unit cell contains 
twelve vertices z/12 = 12, and the coordination number is ki2 = 5. By the choice of the unit 
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shown in 


Fig. 


□ (a) 
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■ (3.1) 



v 

The determinant can be calculated to be 

1)12(^1,^2) = 16{1650732- 680016(cos^i + COS02) + 10151(cos2 ^1 + cos^ ^2) 

-300022 cos 01 cos 02 - (cos^ 0i + cos^ ^2) - 5567 cos 0i cos 02 (cos 0i + cos 02 ) 



+ 158 COS^ 01 COS^ 02 — cos 01 COS 02(cos^ 01 + COS^ 02)} 5 



(3.2) 



such that the numerical evaluation gives 

rf02 



1 rrf0i , 

Zi2 = — / / In 

12 2n 27r 



A2(01,02) 



1.409737903756929. 



(3.3) 



2. Net 13 

Net 13 is obtained by intergrowth of 3^.4.3.4 and 3^.4^ vertices. It can be constructed by 
starting with the square lattice and adding appropriate diagonal edges as shown in Fig. [1] 
(b), such that the coordination number is ki^ = 5. Taking three contiguous squares where 
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two of them contains diagonal edges as a unit cell with 1^13 = 8, we have 
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(3.4) 

The determinant can be calculated to be 

^3(^1,^2) = 4{25560- 2328cos^i - 25872 cos ^2 + cos^ ^1 + 5397 cos^ ^2 

-2158 cos 01 cos 02 - 288 cos^ 02 - 312 cos 0i cos^ ^2 + 4 cos^ 02 
-4cos0icos^e2} , (3.5) 



such that 



1 r d0i r d02 , 



1.409133286424679.. 



(3.6) 



Dvi{0U02) 

The z values of the variants of net 13 shown in Fig. 15 of pLj] are close to this, and are not 
given here to save space. 



3. Net 14 

Net 14 is a simple combination of hexagons and triangles as shown in Fig. [T] (c), such 
that the coordination number is ki4 = 4. Taking a triangle as a unit cell with z/14 = 3, we 
have 



MM, 02) 



4 -1 - e'^^ -1 - e'^^ ^ 



The determinant can be calculated to be 



v 



(3.7) 



Du{0i,02) = 4(9 -3 cos ^1 - 5 cos ^2 - cos ^1 cos ^2) 
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(3.8) 



such that 



1 r dOi r d92 , 

In 



-TT 27r 2tt 



In 4 1 r dOi 



hi 



n 



9 - Scos^i + V8(l - cos^i)(7- cos^i 



In 2 2 r ddi 

h - 

3 3 



— ln(2sin^ + W6 + 2sin2^) . 



2 



(3.9) 



TT V 2: V 2 

An exact closed-form expression for this integral can be derived as follows. After changing 
the variable 9i = 26, we have 

4 r^/"^ 



4 ( /3 1 

Zi4 = In 2 + — y c/^ In sin ^ + y - + - sin^ d 

3 V2' 2) ' 



1 



where 

/(a, 6) 

TT 

and we consider < 6 < 1 < a. For 6 = 0, it can be shown that 

1 /-TT/a 



^ d6' In {^Yn 6 + Va + b sin^ 6* j 



J(a,0) = - / de ln(sin^ + v^) 



1 



-- ln[2(v^ + V^^)] + -Ti2( + V^^) 



TT 



where Ti2(x) is the inverse tangent integral 20 1 

tan-i t 



Ti2(x) 







t 



dt = j2i-^)''' 



X 



2k-l 



k=l 



{2k -ly 



2i 



[Li2(ia;) - Li2(-zx)] 



Here dilogarithm Li2(2;) is defined by 



00 L. 



fc=l 



ln(l -t) 
t 



dt . 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



Taking the derivative with respect to 6 then performing the integral over Q in Eq. (13.111) . 

we get 

d 



db^^""'^^' 2n{l-b)\^ 



fitan"^ 



a + b-1 



tan y/a + b — 1 



1 



4(1-6) 



. (3.15) 
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It can be integrated to give 



I{a,b) = l{a,0) + / r{a,x)dx 



2 



-Tiaf^a + Va - 1) + - ftanh"^ Vb) ftan"^ ^ 



+ — ^tanh ""^ Y 1 tanh ""^ y 1 j tan ^ 



+ 



k=l 



In 

1 

71 V 



y/g - 1 - ^ ^fc ^ ^ y/g + 6 - 1 - ^ ^fc ^ ^ 1 - y/b ^ 
^/a — 1 + y/a + b — 1 + y/g 1 + \/6 

(l-6)(l + v^) 



sin(/c0) 
A;2 



4(2g - 1 + 2y/a{a-l)){l - Vb) 



n 



Ti2(v^ + Va- 1) 



1 - b 



(tanh tanh ""^ \/l tanh ^ tan ""^ i/g 



. Vg — 1 — y/a 



Li2( 



47ri L v'g - 1 + ^g 



e^'^)-Li2( 



-Li2(- 



Vg + &— T + i/g 



+Li2 ^e- 



Li2(- 



v^g + 6 — 1 + -\/g 
-Li2(e'^) + Li2(e- 



(3.16) 



1 + Vb 1 + Vb 

where (p = tan~^ [2^/^/(1 — a)] and tt/2 < (p < n. When b is set to zero, it is clear that 
Eq. ( I3.16P reduces to Eq. (13.12^ . We notice that when a is set to one, the expression for 
J(g = 1, b) can be simphfied, using the identity that Ti2(l) is equal to the Catalan constant 
C = E^=o(-l)"(2^ + 1)"^ = 0.9159655941772190..., as 

1 . /1-V^N 



/(M,4.n(V 



—— ) + - + 7; tanh"^ Vb + -Ti2 - 

V 4 / TT 2 TT VI 



+ Vb 



(3.17) 



which is equivalent to Eq. (28) of [13]. Evaluating I{a,b) in Eq. (13.16^ at g = 3/2, b = 1/2 
and substituting into Eq. (13.101) . we obtain the exact closed-form expression 



1 f3 + 2^/2\ 8^. /V3 + l\ 4 / 



V2 + ^/3/ 37r 



Zi4 = - In , ^ / I 

3 \2 + V3^ 371 



V2 



— ^tan ^ \l —){ tanh 



_i 4^3 -v^- 7 \ 
v^- 2^6 + 7^2/ 



37ri 



U2{iV3 - 2)e^*°) - Li2((y3 - 2)6"^'^°) - Li2((2y6 - 5)e^^°) 
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+Li2((2y6 - 5)e-'*°) + Li2((3 - 2v/2)e''^°) - Li2((3 - 2^pl)e-''^'^) - U^ie''^'') 



+Li2(e^ 



= 1.127778363805542... 
where 0o = tan'i(-2V6) = 1.772154247585227. 



(3.18) 



^. Net 15 

Similar to the Kagome lattice, (3.6.3.6), net 15 is a simple combination of hexagons and 
triangles as shown in Fig. [1] (d). A primitive unit cell contains eleven vertices vi^ = 11. 
There are two 6^, four 3^.6^ and five 3.6.3.6 vertices in each unit cell, so that the effective 
coordination number is = 42/11. Referring to Fig. [1] (d), if one takes the direction to 
the right as 9i and the direction to the upper-right as 62, then 
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(3.19) 



with determinant 



= 4{39910 - 16905(cos^i + cos( 



-^2) 



62) - 1270cos(^i 
5352 cos 9i cos 62 — 25(cos 9i + cos 62) cos(6'i — 62 



- 356(cos=^ 9i 



COS^ 62 



-70 cos 6*1 cos 6*2 (cos 6*1 + cos6'2)} 
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(3.20) 



However, if one still takes the direction to the right as 9i but the direction to the upper-left 
as 6*2, then 
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(3.21) 



with determinant 

^15(^1, O2) = 4{40266 - 16975 cos^i - 1270 cos^s - 16905 cos(^i + 62) - 356 cos^ 62 
-25 cos 61 cos 62 - 5352 cos 6*1 cos(6'i + 62) - 25 cos 62 cos(6'i + 62) 
+70 cos^ 61 + 70 cos 61 cos^ 62 - 70 cos^ 61 cos(^i + 62) 

+712 cos6'i cos 02 cos(^i + ^2) - MOcos^^i cos^^2 cos(6'i + 6'2)} , (3.22) 

which looks distinct from Di^i^Oi, 62). Nevertheless, both determinants give the same asymp- 
totic growth constant and we have the integral identity 

1 r dOi r (162 
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11 2n ./^ 271 



DM,e2) 



^15(^1, ^2) 



1.073270254423056... 



(3.23) 



In contrast to the integral identities in 



141 ] which was obtained by different choice of unit 



cells, here we obtain an integral identity by different choice of directions in the calculation. 
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5. Net 16 



Net 16 is a combination of triangles, squares and hexagons. There are two kinds of net 16, 
namely net 16(a) and net 16(b). Their unit cells are shown in Fig. [I](e) and (f), respectively, 
and both of them contain twelve vertices i^i6(a) = ^^16(6) = 12. For net 16(a), there are six 
3^.4.3.4 and six 3.4.6.4 vertices in each unit cell, so that the effective coordination number 
is /«i6{a) =9/2. We have 
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(3.24) 



The determinant can be calculated to be 

DiQ(a) {Oi, 02) = 112{49746 - 15900(cos 9i + cos 62) - 15899 cos(ei + 62) - 749 cos 9i cos 62 
-749(cos 6*1 + cos 62) cos(6'i + 6*2) + (cos^ 61 + cos^ 62) 



-(cos^ 61 + cos^ 62) cos(6'i + 62) + 204 cos 61 cos 62 cos(6'i + 62) 
-2 cos6'i cos 6*2 (cos 6'i + cos6'2) cos(^^i + 6*2)} , 



(3.25) 



such that 



12 / . 27 
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/^16{a)(^l,^2) 



1.280287248642483. 



(3.26) 
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For net 16(b), there are six 3^.4^ and six 3.4.6.4 vertices in each unit cell, so that the 
effective coordination number is also Ki6(fe) = 9/2. We have 
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v 

The determinant can be calculated to be 

^16(6) (^1,^2) = 16{337746- 109263(cos^i + cos ^2)- 109262 cos(ei-e2) 

-4024 cos 9i cos 62 - 4023(cos 9i + cos ^2) cos(0i - ^2) + (cos^ Oi + cos^ ^2) 

+2118 cos 6*1 cos 6*2 cos(^^i — 62) — sin 61 sin 6*2 (cos^ Oi + cos^ O2) 

—3 cos6'i cos 62 cos(6'i — 6*2) (cos 61 + cos O2) — 2 cos^ Oi cos^ O2 cos(6'i — 62)} , 

(3.28) 



such that 



1 r dOi r (102 , 

which is slightly less than Zi6(a)- 



-Di6(b)(6'i, 62) 



1.277617926708331... 



(3.29) 



6. Net 17 



Net 17 is obtained by intergrowth of 3^.4.3.4 and 3.4.6.4 vertices. A primitive unit cell 
containing twenty-eight vertices 1^17 = 28 is shown in Fig. [1] (g). There are sixteen 3^.4.3.4 
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vertices and twelve 3.4.6.4 in each unit cell, so that the effective coordination number is 
Ki7 = 32/7. Let us write Mir{9i,92) = - M"^{9i,92), where is a diagonal matrix 
with diagonal elements {5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 4, 4, 4, 4, 4, 4, 5, 5, 4, 4, 4, 4, 4, 4, 5, 5} and 
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(3.30) 

The determinant can be calculated to be 

^17(^1, O2) = 16{466754878482464 - 193533953205944(cosei + cos 62) 
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+3411923764300(cos2 9i + cos^^s) - 84153824536936 cos cos02 



-1683339412(cos^ 9i + cos^ 62) - 1185500641948 cos^i cos ^2(cos Oi + cos 62) 
+38809(cos^ 61 + cos ^2) 729383812 COS ^1 COS (cos + cos 62) 
+ 18854974886 cos^ 9^ cos^ 62 - 11032 cos^i cos02(cos=^ 6^ + cos^ 62) 
-11681168 cos^ 9^ cos ^ofcos^i + cos6'2) 

+784 cos^ 01 cos 6'2(cos + COS ^^2) ^ 1568 cos^i cos^2} ; (3.31) 

such that 

= 1.299177753544099... (3.32) 
B. Nets with pentagons 

In this subsection, we consider important nets which involve pentagons. Most of them 
also contain triangles and squares, in addition to pentagons. 

1. Net 18 

Net 18 is a pentagon-only net, and is the dual of the (3^.4.3.4) lattice. The normal 
appearance that it contains equal-sided (but not regular) pentagons with two angles of 7r/2 
is not crucial for the calculation of spanning trees. Let us draw it as shown in Fig. [2] (a) 
with six vertices in a unit cell z/ig = 6. There are two 5^ and four 5^ vertices in each unit 



1 r de, r d92 r i 

28/. 27/^27 ^"[^-('^''^)J 
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cell, so that the effective coordination number is nis = 10/3. We have 
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3 ) 




be calculated to be 










— 36 (cos 9i + cos 62) + 


(cos^ 01 


+ cos^ 


^2)- 


14 cos 6^1 cos 62} 



(3.33) 



which is the same as that for the (3^.4.3.4) lattice as expected. According to Eq. (12.51) for 

= 5, we get 

2 



^18 - -^(3M.3.4) g 



1 r r rf^2 , 

m 

27r 



/^18(^l,^^2) 



0.9405704304962232... (3.35) 



where 2;(32.4,3.4) is given in y^]. A closed-form expression for the integral in Eq. (I3.35P is 
given in [21j] such that 



In 6 AC 

Zl8 = -TT + — + 



371 37T 



Ti2(pg,g) + Ti2(pg, -g) + Tisf-, - 

\q q 



\q q 



-Tia ( pq, -) - Tia (pq, - Tia (-, —) - Tig (-, — '■ 

py V py \q pqy \q pqy i 



(3.36) 



where p = y/S — \/2, q = \f2 — 1, and 'Y\2{x^y) is the generalized inverse tangent integral 



^ tan-i t 
t + y 



dt . 



(3.37) 
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2. Net 19 



Net 19 is shown in Fig. [2] (b), where a primitive unit cell contains nine vertices z/ig = 9, 
and the coordination number is ^ig = 4. By the vertex labeling given in Fig. [2] (b), we have 
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(3.38) 



The determinant can be calculated to be 



ii,02) = 80{442- 191(cosei + cose2)+5(cos2^i + cos 62) — 68 cos 6*1 cos 6' 
— cos 9i cos 62(008 9i + 00862)} , 



such that 



=1 r ^ r 



Dig{6,,62) 



1.144188002944693. 



(3.39) 



(3.40) 
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3. Net 20 



Net 20 is shown in Fig. [2] (c), where a primitive unit cell contains ten vertices V20 = 10, 
and the coordination number is K20 = 4. We have 
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The determinant can be calculated to be 



= 16{7375 - 2995(cos0i + cos^a) + 34(cos2 



(3.41) 



cos^ 62) — 1393 cos 9i cos 62 



— 30 cos6'i cos^2(cos0i + cos^a)} , 



(3.42) 



such that 



'''-To 



2^ 



(162^ 
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27r 



1.150677474300389.. 



(3.43) 



4. Net 21 



A primitive unit cell for net 21 containing twenty vertices z/21 = 20 is shown in Fig. [2] 
(d), and the coordination number is ^21 = 4. Let us write M2i{0i,92) = — M2i{9i,92), 
where M21 is a diagonal matrix with all diagonal elements equal to four and 
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(3.44) 



The determinant can be calculated to be 

D2i{0i, 62) = 16{816275712 - 342967936(cos^i + cos ^2) + 6766760(cos2 61 + cos^ 62) 

- 139108816 cos ^1 cos ^2 - 7200(cos^ + cos^ 62) 

-2389920 cos 61 cos ^2(cos 61 + cos 62) + (cos^ Oi + cos^ 62) 

-1160 cos 9i cos e2(cos2 ^1 + cos^ ^2) + 32014 cos^ 61 cos^ 62} , (3.45) 



such that 



Z21 



]_ r d9i r 
20 2^ 



de2 

2tt 



In 



D2i{ei,e2) 
22 



1.155959257782222.. 



(3.46) 



5. Net 22 



Net 22 contains only triangles and pentagons. Let us draw it as shown in Fig. [2] (e) with 
five vertices in a unit cell 1^22 = 5. There are two 5'^ and three 5^.3 vertices in each unit cell, 
so that the effective coordination number is K22 = 18/5. We have 
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-1 3 
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_g-iei Q 

_1 _e-*(ei+92) 
4 -1 



V -1 

The determinant can be calculated to be 

^22(^1,^2) = 2{99-31(cos^i + cos^2 + cos(ei + ^2)) 

-2(cos^ 



5 6*1 cos 6*2 + (cos 6*1 + cos 6*2) cos(6'i + ^2))} , 



such that 
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1.024172110372259.. 



(3.47) 



(3.48) 



(3.49) 



6. Net 23 

A primitive unit cell for net 23 containing nine vertices z/23 = 9 is shown in Fig. [2] (f), 
and the coordination number is k,23 = 4. We have 
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(3.50) 
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The determinant can be calculated to be 



^23(^1,^2) = 8{4714- 1473(cos0i + cos02 + cos(0i + ^2)) 

-103(cos 61 cos 62 + (cos 61 + cos 62) cos(6'i + 6*2)) 
+ 14 cos 9i cos 62 cos (6'i + 6^2 ) } , 



such that 
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— TV ^ — TT 
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1.152329841150063.. 



(3.51) 



(3.52) 



7. Net 24 

The coordination number for net 24 is K24 = 4. The unit cell given in Fi 
contains ten vertices 1^24 = 10. Using the vertex labeling given in the left-hand 
[2] (g) , we have 
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-side of Fig. 



M24(^l,^2) 



v 



4 


-1 


1 








a 

13 


—a 











-1 


4 


-1 











—a 





-1 







-1 


4 


-1 





—a 




















-1 


4 


-1 











-/3 


-1 











-1 


4 


-1 










-1 


a 





_i 

a 





-1 


4 


-1 











_1 

a 


_1 

a 











-1 


4 


-1 




















1 





-1 


4 


-1 


-1 





-1 





1 

/3 











-1 


4 


-1 











-1 


-1 








-1 


-1 


4 



(3.53) 



with determinant 

£'24(^1,^2) = 4{30135 - 5929 cos ^1-23101 cos ^2 + 36 cos2^i + 2456 cos^ ^2 

-3341 cos 61 cos ^2-4 cos^ 62 - 248 cos 61 cos^ O2 - A cos 61 cos^ ^2} 

/ 2 ^1 2 ^2 ^1 ^2 ^1 s ^2 \ 

= 16(121 — 6 cos -- — 52 cos -- — 59 cos -—cos-- — 4 cos -—cos -— 1 

\ Z Zj Zj Z Z / 



24 



121 — 6 COS — — 52 COS 1- 59 cos — cos 1-4 cos — cos 



2 . A ^1 3 ^2' 



, , , .(3.54) 

2 2 2 2 2 2/^^ 

However, if one takes the primitive unit cell shown in the right-hand-side of Fig. [2] (g) 
containing five vertices, 1^24 = 5, we have 
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(3.55) 



with determinant 

^24(^1,^2) = 2{184-61cos^i -46cos^2 -61cos(^i + ^2) -2cos^iCos^2 



-12cos0icos(^i + ^2) - 2cose2Cos(ei + ^2)} ■ (3.56) 

Although the two determinants given in Eqs. ( 13.54^ and ( 13.56^ are distinct, they give the 
same asymptotic growth constant for net 24, 

-024(^1, 62) 
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10 ./ . 2tx ./ ^ 27r 
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In 



D2i{0ij O2) 



1.148658687301195... 



(3.57) 



5 2n 27r 

This is an example of the integral identity one can obtain by choosing different unit cells in 



the calculation 



8. Net 25 



Let us draw net 25 as shown in Fig. [2] (h) with six vertices in a unit cell z/25 = 6, and the 
coordination number is K25 = 4. We have 
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The determinant can be calculated to be 

£'25(^1,^2) = 4{295 - 104cosei - 136 cos^a + cos^ ^1 + 2 cos^ ^2 - 56 cos^i cos^s 



—2 cos6'i cos^ 62} 



(3.59) 



such that 
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1.150037457106072. 



(3.60) 



C. Nets with heptagons, enneagons or octagons 



In this subsection ,we consider important nets which involve heptagons, enneagons or 
octagons. 



1. Net 26 

Net 26 contains only squares and heptagons. A primitive unit cell containing twelve 
vertices 1^26 = 12 is shown in Fig. [3] (a), and the coordination number is K26 = 3. We have 
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The determinant can be calculated to be 

^6(^1, ^2) = 4{4140 - 1752(cos0i + cos^2) + 37(cos2 0i + cos2 02) - 686cos^iCos^2 
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- 12 COS ^1 008^2(008^1 +008^2)} , (3.62) 



8uch that 




(c) (d) (e) 

FIG. 3: (a) A unit cell of net 26. (b) Net 27. (c) A unit cell of net 28. (d) A unit cell of net 29. (e) Two 
unit cells of net 30. Vertices within a unit cell are labeled. 
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2. Net 27 



Net 27 contains only triangles and enneagons as shown in Fig. [3] (b), where each unit cell 
contains four vertices z/27 = 4, and the coordination number is K27 = 3. We have 

/ 3 -1 -e'^^ -e-*^i \ 
-1 3 -1 -1 
-e-'^^ -1 3 -1 

3 / 



\ -e*^i -1 -1 
The determinant can be calculated to be 



(3.64) 



^27(^1,^2) = 8{3-(cos^i + cos^2)-cos(^i + ^2)} , 



which is four times of that for the triangular lattice 
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0.7504060243042857... 



(3.66) 



3. Net 28 



Net 28 is the B net of YCrB4, which contains only pentagons and heptagons. A primitive 
unit cell containing sixteen vertices z/28 = 16 is shown in Fig. [3] (c), and the coordination 
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number is K28 = 3. We have 
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(3.67) 

The determinant can be calculated to be 

^23(^1,^2) = 4{ 12 1795 - 115860 cos 01 - 11428 cos ^2 + 19777 cos^^i + 2 008^62 



- 1 1472 cos 9i cos ^2 - 628 cos^ - 2 1 32 cos^ 9i cos 02-2 cos 9i cos^ 62 



+4 cos^ 9i - 56 cos 6'icos6'2} 



(3.68) 



such that 
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0.8025881723106822. 



(3.69) 



4. Net 29 



Net 29 is the B net of Y2LnB6, which contains pentagons, hexagons and heptagons. A 
primitive unit cell containing twenty four vertices V29 = 24 is shown in Fig. [3] (d), and the 

29 



coordination number is K29 — 3. Let us write M29{9i, 62) — - 
a diagonal matrix with all diagonal elements equal to three and 
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M2g{0i,92), where M^g is 
\ 



V 



The determinant can be calculated to be 



^29(^1,^2) = 8{36407767 -21392474 cos ^1 - 9913039 cos ^2 + 1035241 cos^^i 

+47642 cos^ 02 - 5861590 cos 9i cos 62 - 4068 cos^ 9i - 288591 cos^ 9i cos 62 
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-31198 cos 01 cos^ 02-2 cos^ 02 + 2 cos^ 0i - 840 cos^ 0i cos 02 



such that 



+1156 cos^ 01 cos^ 02-2 cos 0i cos"* 02-4 cos"* 0i cos^ 02} 



= ^ r ^ r ^ In 



-029(01, 02) 



0.8043880179770491. 



(3.71) 



(3.72) 



5. Net 30 

Net 30 is a net with 5^.8 and 5.8^ vertices. Its coordination number is K30 = 3. The unit 
cell given in Fig. 31 of [1)] contains twelve vertices v^q = 12. Using the vertex labeling given 
in the left-hand-side of Fig. [3] (e), we have 



M3o(01,02 



3 


-1 








-1 


_1_ 

a 

















\ 


-1 


3 


-1 


a 





























-1 


3 


-1 




























—a 


-1 


3 


-1 























-1 








-1 


3 


-1 




















—a 











-1 


3 


-1 
































-1 


3 


-1 











—a 




















-1 


3 


-1 








-1 























-1 


3 


-1 


—a 











1 

/3 

















-1 


3 


-1 





























_i 

a 


-1 


3 


-1 




















_1 


-1 








-1 


3 > 



(3.73) 



v 

with determinant 

D-io{9i, 02) = 4{5160 - 6791 cos 0i - 169 cos 02 + 2261 cos^ 0i - 221 cos 0i cos 02 
-184 cos^ 01 - 56 cos^ 0i cos 02 + 4 cos^ 0i - 4 cos^ 0i cos 02} 

/ 2 ^1 01 02 4 01 ■? 01 02 

= 16 60 — 49 COS — + 11 cos — cos — + 4 cos — +4cos — cos — 



31 



60 - 49 cos^ ^ - 11 cos ^ cos ^ + 4 cos^ ^ - 4 cos^ ^ cos ^) . (3.74) 
2 2 2 2 2 2/^^ 

However, if one takes the primitive unit cell shown in the right-hand-side of Fig. [3] (e), where 
each unit cell contains six vertices, z/30 = 6, we have 
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(3.75) 



with determinant 

^30(^1, O2) = 2{73 - 45 cos 9i - 13 cos 02 - 13 cos(^i - 62) + 2 cos^ 9i - 2 cos 9i cos 62 



-2cos^iCos(^i -62)} . 



(3.76) 



Although the two determinants given in Eqs. (13.741) and (I3.76P are distinct, they give the 
same asymptotic growth constant for net 30, 







12 













2tt 



In 



D 



2n 



In 
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0.7985013545791521... 



(3.77) 



This is another example of the integral identity one can obtain by choosing different unit 
cells in the calculation jl^ . 



IV. DISCUSSION 



It is of interest to see how close the exact results presented above are to the upper bound 
given in Eq. (12. 8p . For this purpose, we define the ratio 



In bk 



(4.1) 



32 



for a /c-regular lattice Ak, where bk is given by Eq. (12. 7p . For a lattice A which is not 
fc-regular, we replace /c by k in Eq. (12.71) and consider the ratio 



The values of and for various 



TA = ^ . (4.2) 
In 6k 

attices A are summarized in Table [H Our results 



agree with the observations made in 1^ that za is relatively large for large value of k (or 
k). For the two-dimensional lattices with k = 4 studied here and in 14 1, their values 



of za are all smaller than that of the square lattice, Zsq = AC/n = 1.166243616123275 
which indicates that the square lattice may be the most densely connected two-dimensional 
lattice with k = 4 [l^ . Similarly, the values of z\ for the lattices with k = 3 studied here 
and in 



12 



Il4| are smaller than that of the honeycomb lattice, Zhc = 0.8076648680486262... 
{s], which indicates that the honeycomb lattice may be the most densely connected two- 
dimensional lattice with k = 3. The triangular lattice, the dual of the honeycomb lattice, 
would be the most densely connected two-dimensional lattice with k = 6. 

The ratios r\ are close to each other no matter the lattice A is fc-regular or not. We 
thereby conjecture the upper bound 

ZA < ln(6,) (4.3) 

for a lattice A with effective coordination number k > 3, which generalizes the bound given 
in Eq. (12. 8p for /c-regular lattices. 
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